Abstract. We present ISO-LWS measurements of the 
main gas cooling lines, [0+] 158/Ltm and [0°] 63/^m towards 
a moderate opacity molecular cloud {Av ~ 3), L1721, illu- 
minated by the B2 star 1/ Sco (x = 5-10). These data are 
combined with an extinction map and IRAS dust emis- 
sion images to test our understanding of gas heating and 
cooling in photo-dissociation regions (PDRs). This nearby 
PDR is spatially resolved in the IRAS images; variations 
in the IRAS colors across the cloud indicate an enhanced 
abundance of small dust grains within the PDR. A spa- 
] tial correlation between the gas cooling lines and the in- 
■ frared emission from small dust grains illustrates the dom- 
' inant role of small dust grains in the gas heating through 
• the photoelectric effect. The photoelectric efficiency, de- 
O ^' termined from the observations by ratioing the power ra- 
diated by gas and small dust grains, is in the range 2 to 
3 % in close agreement with recent theoretical estimates 



( [Bakes fc Tielens, 1994| ; [Weingartner fc Draine, 2001| ) 



< 

I I The brightness profiles across the PDR in the [C+] 158/im 
^ . and [O*'] 63/im lines are compared with model calculations 
' where the density profile is constrained by the extinction 
data and where the gas chemical and thermal balances are 
solved at each position. We show that abundance varia- 
tions of small dust grains across the PDR must be consid- 
ered to account for the LWS observations. 
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1. Introduction 

The bulk of interstellar matter is found in regions of low to 
moderate opacity to UV and visible light where stellar ra- 
diation plays a dominant role in determining the chemical 
and thermal state of the gas. These photon-dominated or 
photo-dissociation regions have been the subject of many 
observations and theoretical studies over the past twenty 
years (Hollcnbach & Tielens, 1999, for a review see ). In 



calculations give a net estimate of the [C+] 158/im line 
emission per hydrogen atom in agreement with the value 
derived from the comparison of the FIRAS and H I data at 
high Galactic latitude (Boulanger et al., 1996a). However, 



the classical theoretical paradigm, a PDR is character- 
ized by its proton density, uh, and a scaling factor, Xj 
which normalises the incident radiation field to the Solar 
Neighbourhood radiation field in the far- ultraviolet. Com- 
parison between observations and model calculations have 
concentrated on bright objects with high x > lO'^ and gas 
density nn > 10^ cm~^ which are the easiest to observe. 
The gain in sensitivity provided by the Infrared Space 
Observatory (ISO) for gas line observations now permits 
to extend these studies to less excited PDRs with lower 
uh and x which are the most widespread in the inter- 
stellar medium (ISM). Models show that in such PDRs 
the gas heating is dominated by the photoelectric effect 
on d ust grains (prlollenbach et al., 1991 : Lc Bourlot et al 



this validation of the BT work depends on the relative gas 
mass in the cold and warm phases of the ISM which is 
poorly constrained observationally. 

In this paper, we investigate the photoelectric heating 
of the ISM through the study of a specific, low-excitation 
PDR. The comparison of the power radiated by gas and 
dust allows us to correlate the gas and dust emission and 
determine the photoelectric efficiency. 

To carry out this study, we selected the Lynds Dark 



Nebula 1721 (Lynds, 1962, L1721, ), a nearby molecular 
cloud in the p Ophiuchi region. This cloud as seen in the 
IRAS data has a roughly spherical geometry and a mod- 
erate opacity to stellar light {Av = 3 at cloud center). It 
is heated on one side by the B2 star v Sco and by a more 
isotropic interstellar radiation field created by the Upper 
Scorpius association. We present in this paper observa- 
tions of the [C^] 158/xm and [0°] 63/im lines obtained with 



the Long Wavelength Spectrometer (Clegg et al., 1996, 
LWS, ) on board the Infrared Space Observatory ( pK^esslei 



the subject uf lecenL llieuicLical iuvesLigaLiuus which piu 
vide an estimate of the photoelectric yield as a function 
of grain size where the small grains (radius <100 A) dom- 
inate the overall photoelectric heating (Bakes & Tielens 



^;|Kcmpci cL al., 1999|). PliutuelcLLiit lieaLiug has been ^^^a j^q ^ ^^^^^ observations are combin^d^ 



1994; Weingartner & Draine, 2001). It is thus in low exci 



tation PDRs including the diffuse interstellar medium that 
one can most directly validate this theoretical understand- 
ing of a key physical process for the ISM, independently 



of other heating processes. Recently, Wolfirc et al. (1995) 
have used the work of Bakes and Tielens (1994, hereafter 
BT) to compute the equilibrium thermal states of the dif- 
fuse atomic medium. For the Solar Neighbourhood, their 
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IRAS images and an extinction map (Sects. 2-3). Thanks 
to its proximity to the Sun and the moderate gas den- 
sity, the L1721 PDR is spatially resolved by both the gas 
and dust observations. The data is interpreted within the 
framework of a model of the gas emission (Sect. 3). The 
data allows us to determine the photoelectric efficiency 
and to illustrate the dominant role of small dust grains in 
the photoelectic heating (Sect. 4). The paper conclusions 
are summarised in Sect. 5. 

2. Observations 

The L1721 cloud is centered at a = 16h 14mn 30.41s and (3 
= -18° 58' 7.9" (Epoch 2000), at the north of the p Ophi- 
uchi complex. The B2 star i/ Sco is located at a distance 
of ^^134 pc. The projected distance between the center of 
the cloud and the star is approximately 1.7 pc. We discuss 
in detail the incident radiation field of L1721 in Sect. 3.2. 
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Fig. 1. Ln21 as seen in the IRAS data. Top: map of the 
sum of IRAS bands (12+25+60+100 fim) with the beams 
of LWS observations projected on sky. Bottom: map of 
the 12/100 fim IRAS band ratio. In this latter map, the 
limb brightening or halo of the cloud is clearly visible. The 
circles show the LWS beam in the cut taken across the 
cloud ( 20 positions ) and the position of v Sco is marked 
by a star. 



of glitches, all scans are again co-added to produce the fi- 
nal line spectrum. Line fluxes and error bars are obtained 
by fitting a Gaussian and a linear baseline to the profiles. 
The Fig. || shows the line intensity profiles. 

To measure the line emission from LI 721, we need to 
subtract the contribution of the diffuse ISM along the line 
of sight. We express the background contribution to the 
line emission as: 

Ib = Yx yl, (1) 
where Y is the line-to-continuum ratio and vli, the con- 
tinuum brightness at the line wavelength. In the case of 
the line, we take the line-to-continuum ratio measured 
in the Galaxy at high latitude, Yc+ = 0.009±0.0015 from 



Boulangcr ct al. (1996a). This value agrees well with the 
Galactic plane mean value of 0.0082 reported by Wright 



et al. (1991) based on the FIRAS (Far-Infrared Abso- 



On the map of the 12/100 jim. IRAS band ratio (see 
Fig. |l|), the limb brightening of the cloud is clearly visi- 
ble. Former studies have shown that ratios of the differ- 
ent IRAS bands mostly trace dust abundance variations 
( Boulanger et al., 1990| ; Bernard et al., 1993 ). We discuss 
the abundance variations of small grains across LI 721 in 
Sect. 3.5. 

2.1. Gas cooling lines 




Fig. 2. Intensities of the and 0'^ lines with 3cr er- 
ror bars as measured by LWS at each position across 
L1721. The background contribution has been subtracted 
(see text). The v Sco star lies at the right end of the cut. 

The [O"] 63/im and [C+] 158/im lines have been ob- 
served by ISO with the LWS (TDT = 45401125, LWS02 
grating mode, Clegg et al., 1996). At each sky position, 
the C"*" and 0° lines have been observed 8 and 30 times 
respectively with 5 seconds-individual scans (the integra- 
tion per grating point was 0.4 seconds). The resolution 
(FWHM) is AA = 0.5 and 0.29 ^lul for the C+ and 0° 
lines respectively and each resolution element was sam- 
pled twice. We use the LWS pipeline version 7 for standard 
data reduction. To remove glitches (cosmic ray impacts on 
the detectors), we co-add all scans to obtain a first esti- 
mate of the mean line spectrum. Subtracting this mean 
spectrum from each individual scan, we obtain a (noise -|- 
glitch) - spectrum in which glitches are identified as 3cr- 
and 5(T-outliers for C"*" and respectively. After rejection 



lute Spectrophotometer) survey of the far-infrared emis- 
sion from our Galaxy. [Wright ct al. (1991 ) observed, as 
expected theoretically, that the intensity of the 0+ line 
scales linearly with the dust continuum emission. The vari- 
ations of the line-to-continuum ratio, due to differences 
in the dust properties, in the color of the ambient inter- 
stellar radiation field or the saturation of the C+ cooling 
rate when the level populations are in equilibrium with 
dense warm gas, are within 25% in the plane of Galaxy. 
In the following we consider that 1^+ — 9±2.25 10^'^ in 
the proximity of the L1721 cloud. For the 0° line, the 
combination of the observations of Caux fc Gry (1997 ) in 
the diffuse ISM and the IRAS band at 60 /xm provide a 
large range of Yqo values from 5 10~* to 5 10~^. To derive 
the 0° background around the L1721 cloud we take the 
median of these values Yqo = 1.7 10""^. Note however that 
the large uncertainty on Yqo will not affect our estimate of 
the O" line emission in the L1721 cloud because the con- 
tribution of the diffuse ISM to the 0° emission along the 
line of sight is much lower than the line intensity in the 
cloud. For the continuum brightness, we use the DIRBE 
(Diffuse Infrared Background Experiment) data corrected 
for zodiacal light. The continuum values are found from 
the average of DIRBE points in the diffuse ISM around the 
L1721 cloud. The variations of the DIRBE values are typi- 
cally within 20% in the surrounding of the cloud. We find: 
z^^(158 Atm) = 9.2±1.9 10"'' ergs/s/cmVsr and vl^{m 
/^m) = 3.4±0.7 10~^ ergs/s/cm^/sr. The corresponding 
background emission values which have been subtracted 
from our LWS cut are 8.3±2.7 10"^ and 5.8±1.2 lO"'^ 
ergs/s/cm^/sr for C+ and 0° respectively. 

2.2. Dust infrared emission 

Fig. H shows the background subtracted surface brightness 
in the four IRAS bands. For each IRAS band, the back- 
ground emission has been estimated along cuts extending 
the LWS cut at both ends of the cloud. The background 
emission is taken to be the average of the minimum values 
of these extended cuts: = 1.2±0.2, /fg = 1.8±0.3, = 



4 
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Fig. 3. Surface brightness in the 4 IRAS bands as a func- 
tion of the observed LWS positions. The background emis- 
sion has been subtracted (see Sect. 2.2). 



7 ± 1.6 and /fpo = 35 ± 3.9 MJy/sr. 

The IRAS and COBE data are well explained with dust 
models comprising three components (Desert et al., 1990; 
Dwek et al., 1997). We adopt the Desert et al. model: the 



three dust components are thus by order of increasing sizes 
(in terms of the grain radius, a): 

1. large aromatic molecules or PAHs (Polycyclic Aro- 
matic Hydrocarbons) which contribute most of the 12 
/zm band (a = 4 to 12 A) 

2. the Very Small Grains (VSGs) which dominate the 25 
fim and 60 /im bands (a — 12 to 150 A) 

3. the Big Grains (BGs) which make most of the 100 fim 
band (a = 1.5 10"^ to 0.11 fim). 

In the following, we will often refer to PAHs and VSGs 
as the small dust grain populations. To derive the photo- 
electric efficiency for each grain component, we need to 
decompose the total infrared (IR) emission into individ- 
ual dust contributions at each position in the cloud. 

In the case of PAHs, we use the foll owing relationship 
based on ISOCAM and DI RBE data ( |Boulanger et al 
1996b| ; [Bernard et al., 19941 ): 



I PAH — 1.5xi'I^{12 /im) erg cm 



sr 



The emission of big grains is modelled with a modi- 
fied blackbody of temperature Tbg and an emissivity law 
proportional to i/^. The value of Tbg is estimated using 
the flux measured by all ten LWS detectors covering the 
range A = 42.5-176 /zm. At each position, we fit the ten 
spectral points with the modified blackbody: a constant 
temperature of Tbg = 20 K at all positions is found to 
provide a satisfactory fit. Finally, this modified blackbody 
is scaled to match the 100 /im IRAS band. 

For the VSGs emission, we first obtain two spectral 
points at 25 and 60 /xm from: 

hsMVSG) = I25.MIRAS) - I25MBG) 

At each cloud position, we fit these two points with 
a modified blackbody of constant temperature TysG = 
70 K and an emissivity law proportional to v. 

The bolometric intensity of each dust component is 
shown in Fig. ^: in fact, the VSGs and BGs emission pro- 



files are quite similar to the 25-1-60 fim and 100 /im IRAS 
band profiles respectively. We note that BGs and VSGs 
temperatures are expected to decrease as one penetrates 
into the cloud. The constant temperature we find may be 
due to the spherical geometry of the cloud where warm 
and cold dust are mixed along a given line of sight. 




Fig. 4. Integrated flux of each dust component as a func- 
tion of the position along the LWS cut. v Sco lies at the 
right end of the cut. 



3. Modelling the gas line cooling 

To study the photoelectric effect on dust grains across this 
nearby PDR, we will compare the total gas cooling lines 
to the dust IR emission (see Sect. 4). In the case of low ex- 
citation PDRs, cooling by the pure rotational lines of H2 
may be important. As this latter has not been observed, we 
resort to an updated version of the PDR model described 
in Le Bourlot et al. (1993) to estimate the cooling by H2. 
Moreover the comparison between the brightness profiles 
observed across the L1721 PDR in the [C+J 158/xm and 
[0°] 63/im lines with model calculations provides a test for 
our understanding of gas heating/cooling in PDRs. The 
key ingredients of this model are the density profile across 
the PDR and the incident radiation field at the cloud sur- 
face. In the following, we describe how these quantities 
are constrained in order to produce a realistic model of 
the L1721 PDR. 



3.1. The density distribution 

To determine the density profile of LI 721 we use a vi- 



sual extinction map derived from star counts (Cambresy, 
|l999|) . Fig. H displays the column density, Nh = 
Av/{0.54: X 10^^^) cm~^ observed across L1721. The col- 
umn density profile is well represented by a gas density 
distribution of the form (see Fig. ||): 



nnir) 



n% X {r/riY 



'■H 



ri < r < r„ 
r < ri 



(2) 



with r the radial distance to the cloud center, n^ = 2915 
cm^'^ the constant density in the < r < ri region with ri 
= 0.2 pc and rm = 2.2 pc. In this derivation, we assume 
that Ry is constant and equal to 3.1 across the L1721 
cloud. A higher value of Ry, 5.5, that is commonly ob- 
served in opaque cores, would multiply our extinction val- 
ues by a small factor, about 1.12. 
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Fig. 5. Comparison of the column density observed (solid 
line) across L1721 and the analytic function fitted to it 
( dotted line ) as a function of radial distance, r, and the 
visual extinction integrated from the cloud border, Ay . The 
peak observed at the right end of the cut is due to the 
plate saturation by the v Sco star. The LWS positions are 
marked by squares. 



3. 2. The radiation field 

The ratio /looMm/^H provides a good mea sure of the ra- 
diation fi^ld^trengthjBema^^^^ 



Fig. 7. Geometry used for modelling L1721. The projected 
distance of v Sco to the cloud center is 1.7 pc. 

on the position of the star along the hne of sight: we refer- 
ence this position with an angle 9 between the sky plane 
and the line connecting the cloud center to v Sco (see Fig. 

Next, the ISRF at this position is likely to be enhanced 
relative to the mean Solar Neighbourhood ISRF because 
L1721 is heated by the B and A stars cluster of Upper 
Scorpius. The Upper Scorpius complex contains a large 
proportion of stars with effective temperature of less than 
20,000 K ( |dc Geus et al., 198S| ). We use the [Mathis et al" 



In Fig. t, we show the relationship between the IRAS 100 (1983) radiation field to represent the mean Upper Scor 



/im surface brightness and the observed column density. 
In the diffuse part of the cloud {Ay < 0.8), we observe a 
good linear correlation between the IRAS 100 fim surface 
brightness and the column density. The slope hoofmi/NH 
of this correlation provides an estimate of the radiation 
field intensity. For a set of clouds in Chamaleon with an 
opacity comparable to L1721 Boulanger et al. (1998D found 



/Nh ratio of 7.5±2.5 MJy/sr/mag which corre- 
sponds to 4±1.3 MJy/sr/(10^^ cm~^). This study provides 
the reference value for the mean Solar Neighbourhood 
ISRF (IntersteUar Standard Radiation Field). In LI 721 
we fin 4 a slope of 15 ±4 for the cloud side opposite to 
Sco and 23 ±4 for the cloud side facing Sco (see Fig. ||) 
which corresponds to 4 and 6 times the ISRF respectively. 
We emphasize here that the IRAS 100 /im emission pro- 
file is relatively insensitive to abundance variations of the 
small grains (PAHs and VSGs) as shown in Bernard et al 
(1992| ). 



plus radiation field. This last radiation field is assumed 
to be isotropic and in the following we call this contribu- 
tion the isotropic radiation field. The radiation field excit- 
ing LI 721 thus only depends on the distance of ly Sco to 
the cloud center and on the strength of the Upper Scor- 
pius contribution. The slope of the 100 fjjm-Nu correlation 
measures the power absorbed by Big Grains per Hydro- 
gen atom. This power depends on the intensity but also on 
the spectral distribution of the radiation field. In practice, 
for a given radiation field intensity at 1000 A, the emis- 



sion of BGs per H atom is 2.5 times larger for the Mathii 
et al. (1983| ) radiation field than for a 20,000 K blackbody. 




Fig. 6. The 100 /im IRAS surface brightness (a) opposite 
to V Sco and (b ) facing v Sco as a function of the column 
density observed along each line of sight. 

The radiation field heating L1721 is the sum of two 
contributions. The first contribution is from the i/ Sco star, 
which has a luminosity of 6300 Lq and an effective tem- 
perature of 20,000 K. The contribution of ly Sco depends 



We have taken into account this dependence to translate 
the 100 fim/Nn slopes into x-values which represent the 
strength of a radiation field at A = 1000 A in units of 
1.6 10""^ ergs/cm^/sr (Habing, 1968). The radiation field 
intensity is described in this way because the far-UV pho- 
tons are more efficient to heat the gas. The strength of 
the v Sco radiation field is represented by Xsco and the 
strength of the isotropic radiation field by the factor xi- 
We have considered two cases: 

1. X/ = 1 for the isotropic radiation field and Xsco — 12.5. 
In that case most of the radiation field strength on 
both cloud sides is due to z/ Sco. 

2. The Upper Scorpius isotropic radiation field on L1721 
is the dominant radiation field on the cloud side oppo- 
site to v Sco and then xi — ^ and Xsco = 5. 

A detailed modelling of the dust IR emission in LI 721, 



based on the model of Bernard et al. (1992), is currently 
underway in our group to discriminate between these two 
cases. We expect that the case 2 radiation field will propa- 
gate more efficiently than the UV light towards the cloud 
center because of the strong visible component of the ra- 
diation field. Also, note that the effects of dumpiness can 



6 



Habart et al.: Photoelectric efficiency across L1721 



affect the penetration of star light inside the cloud. For 
the line modelling, we select the case 2 to describe the 
exciting radiation field of L1721. In Sect. 3.4, we will dis- 
cuss the sensitivity of the gas cooling lines emission to the 
incident radiation field. In the case 2, the contribution of 
V Sco amounts to a mean value of Xsco = 5 for the cloud 
side facing v Sco, which corresponds to 9 = 70° for the 
position of the star. As the far-infrared extinction is small, 
a very similar 100 ^.m emission band profile is expected 
for 9 = -70°. We will see in Sect. 3.4 that 9 = +70° is 
favoured. 



3.3. PDR model 



3.4- Spatial distribution of line emission 

We now present results of the PDR model for the pa- 
rameters described in the last three sections. The total 
strength, xt = Xi + Xsco, of the radiation field on the cloud 
surface ranges from 4 (opposite to z/ Sco) to 11 (facing ly 
Sco). Figs. H & P describe the chemical stratification and 
the thermal budget of the plane-parallel model as a func- 
tion of depth into the cloud. In these figures we have used 
X = Xt = 5, which corresponds to the mean radiation field 
received by a particle located on the cloud surface facing 
I' Sco. We note that the photoelectric effect dominates by 
far the gas heat budget. 



To interpret the gas line emission we make use of the 
model described in Le Bourlot et al. (1993). In this model 



a PDR is represented by a semi-infinite plane-parallel slab 
with an isotropic radiation field incident on the interface. 
The inputs parameters are (i) X; the scaling factor for the 
radiation field, and (ii) the density profile. With these in- 
puts the model solves the chemical and thermal balances 
starting from the slab edge. The scattering of the radia- 



tion by dust grains is treated with the formalism of Flan 



nery et al. (1980) in the plane-parallel case and assumes a 
mostly forward scattering (g =< cosa >=0.6, where a is 
the scattering angle). For extinction properties, we use the 
analytical fit of pitzpatrick fc Massa (1988 ) to the Galac- 



tic average extinction curve. The radiative transfer in the 
absorption lines of H2 and CO is treated in detail and the 
individual line profiles are treated with the prescription of 
Federman et al. (1979 ). We use recent gas phase elemental 
abundances measured in the diffuse interstellar medium: 
He/H = 0.1, C/H ^ 1.4 10'^ (|Cardcm ct al., 1996] ), and 
0/H = 3.19 lO-'' ([Meyer et al., 1998|). 



For the photoelectric effect on small dust grains, we adopt 
the formalism of BT. In particular, our model takes into 
account the actual spectral distribution of the radiation 
field to compute the heating rate (Eqs. 1 and 14 of BT). 
The exponent of the power law size distribution of small 
grains is -3.5 while the lower and upper limits of the grain 
radius are 4 and 100 A respectively. Larger grains are 
not included in the computation of the photoelectric heat- 
ing rate. They have a low photoelectric efficiency due to 
the large number of collisions the electron must undergo 
before escaping the grain (Vcrstraete et al., 199C). The 
computation of the photoelectric heating assumes that the 
small grains are graphitic and spherical. Grains with less 
than 50 C atoms are likely to be planar: these species 
(PAHs) are less charged and consequently have a higher 
photoelectric efficiency (BT). However, fo r a given radi 



at ion field, the gain does not exceed 25% ( Wolfire et al 



1995 ): such variations fall within our error bars on the gas 



cooling lines. 



Fig. 8. Variation of the gas temperature, the gas density 
and the fractional ionization x — n^jnn for a plane- 
parallel PDR model with x — 5 and the density profile of 
Sect. 3.1. We also show the H/H2 and /CP /CO tran- 
sitions for the same model. The v Sco star lies to the left 
end of the figure. 




Fig. 9. Main heating and cooling rates for the PDR model 
(X - 5). 

We now use the PDR model to compute profiles of line 
fluxes across the cloud which can be compared to the ob- 
servations. To do this, we approximate the spherical shape 
of the cloud by a combination of plane-parallel models. 
This allows us to take into account the variation of the 
radiation field at the cloud surface and the fact that each 
line of sight crosses diffuse, warm, interface gas as well 
as more deeply embedded, colder regions. We thus divide 
the cloud into angular sectors. The first sector takes up 
all the cloud side opposite to Sco (the area noted A in 
Fig. In addition, six equal sectors are defined on the 
side facing v Sco (the area noted B in Fig. |^): these sec- 
tors are symmetrically distributed with respect to the line 
connecting v Sco to the cloud center. To each angular sec- 
tor is associated one value of Xt- The three sectors facing 
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V Sco on one side of the v Sco-cloud center line have Xt 
= 11, 10 and 9 (going away from the f Sco-cloud center 
line). The large sector A opposite to v Sco receives Xt — 
Xi = 4. We then run a PDR model for each sector with 
the density distribution of Sect. 3.1 and the corresponding 
value of Xt- The model output is a local line emission rate 
(Jy) which only depends on the optical depth measured 
along a radius from the cloud edge. 




Fig. 10. Geometry used to calculate the line emission pro- 
files intensities. Index k labels the various angular sectors 
while i labels the -layers of the PDR models. 

The line intensity integrated along a given line of sight 
is written as (see Fig. 
^ ^A^{k,i) 



k,i 



An 



X l{k, i) 



(3) 



In this expression, k labels the angular sectors and i the 
various layers (defined by the radial Ay taken from the 
cloud edge) crossed by the line of sight within an angular 
sector; l{k,i) is the length in cm of the layer i in the k 
angular sector on the line of sight. The line cooling rate in 
erg/s/cm'^ A^{k, i) is equal to Jy(fc, i) x (i{i) where Ju{k, i) 
is the line emissivity extracted from the PDR model and 
the escape probability from layer i to the cloud edge 
along the line of sight. For (3, we use the formalism of 
[Tielens fc HoUcnbach (1985| ) in their Appendix B with a 
turbulent Doppler width 5vd — 1 km/s. We find that the 
[C+] 158/im and [O"] 63/xm lines reach an optical depth of 
1 at Ai,=0.5 and 0.3 respectively. The H2 rotational lines 
are optically thin. 

The model line intensity profiles results are compared 
to the data in Fig. 11, We have put the star at = +70° 
(Fig. 0). Actually, for 9 = -70° the model hne emission 
is significantly lower due to line opacity. With a dust- 
to-gas mass ratio of the small grain populations (PAHs 
and VSGs) constant throughout the cloud and equal to 
8.6 10"'', the observed C+ and 0° line intensity pro- 
files are not reproduced by the PDR model. Indeed, the 
[C+] 158/Ltm and [O"] 63/im line intensities are underesti- 
mated by a factor of about 2 and 3 respectively. The un- 
certainty on the radiation field, mentionned in Sect. 3.2, 
amounts to relative error of 40% for C+ and 70% for 0°. 
Such error bars cannot explain the discrepancy between 
the observed and predicted profiles. Also a higher C and 
O abundances cannot accommodated this difference, be- 
cause the fine structure line intensities are practically lin- 




Fig. 11. Line emission profiles as observed (solid lines 
with squares) and predicted by the PDR model as a func- 
tion of radius and of visual extinction from the cloud edge. 
The dotted and dashed lines correspond to a calculation 
with a small dust grain abundance constant throughout the 
cloud and Sv^ equal to 1 km/s and 3 km/s respectively. The 
model with small dust grain abundance variations and Svd 
~ 1 km/s or 3 km/s is represented by the solid and dot- 
dash lines respectively. The H2 line cooling has not been 
observed. 



early dependant on the gas phase elemental abundances. 
Such discrepancies may arise from transfer effects on the 
emerging cooling lines and/or an underestimated heating 
rate. 

In the PDR model, the non-thermal velocity fields are 
represented by the turbulent Doppler width 5vd which is 
added to the thermal width ( Tielens fc HoUenbach, 1985 ). 
This parameter affects the transfer of resonant UV pho- 
tons through the absorption lines of H2 and CO and the 
opacity of the emerging IR cooling lines. We show the 
result of varying Svd in Fig. |ll|: for a reasonable range 
of Jwd-values, we see that decreasing the opacity of the 
emerging cooling lies does not account for the observed 
intensities, in particular [O*^]. 

Alternatively, the low intensity of the gas cooling lines 
in the model may tell us that the gas temperature is too 
low (see Fig. ^). Indeed, the predicted [0°] 63/xm line emis- 
sion is very sensitive to the gas temperature because the 
energy of the transition corresponds to a temperature of 
228 K. Thus, another solution to match the data is to in- 
crease the heating rate. This latter scales with the UV ra- 
diation field intensity (x) , the small grains abundance and 
the photoelectric efficiency (BT). The band ratios shown 
in Fig. |l2| suggest that the PAHs and VSGs abundance 
vary across the cloud. In the next subsection we estimate 
these variations and show that they have an important im- 
pact on the emerging gas cooling. We have not attempted 
to change the abundance of BGs, because as mentioned in 
the Sect. 3.3, these large grains have a low photoelectric 
efficiency. 



3.5. Abundance variations of the small grains 

In this section, we first use the dust brightness profile to 
quantify the abundance variations of small dust grains. We 
then include these abundance variations in the PDR model 
and estimate the corresponding line intensity profiles as in 
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the last section. Here we test the idea that the discrepancy 
between the observed and the predicted profiles is due to 
abundance variations of the small dust grains. Ratios of 
the brightness profiles of Fig. ^ are presented in Fig. 




Fig. 12. Emission ratios of the different dust components 
as observed (solid lines with squares) and as predicted 
( solid lines ) as a function of radius and of visual extinc- 
tion from the cloud edge. 




Fig. 13. Variation of the interstellar carbon abundance in 
PAHs (solid line) and VSGs (dotted line) of PDR models 
as a function of radius and of visual extinction integrated 
from the cloud border. 

We have estimated the abundance variations of small 
grains from the ratio of PAHs and VSGs emission to that 
of BGs. A coherent modelling of the dust emission from 
a spherical cloud (e.g. Bernard et al., 1992, 1993) which 
includes consistent dust abundance variations is beyond 
the scope of the present paper. We estimate the model 
IR emission of dust grains I^Xh^ ^vsg and I^^'^ as fol- 
lows. The power in the IR emission of the different dust 
components, X ~ [PAH, VSG, BG], is equal to the power 
absorbed, i.e.: 

k i 

where 



All 



P. 



13.6 cV 



ax{E)Fk4E)dE 



(5) 



is the power absorbed in erg s^^ H^^ with ax, the ab- 
sorption cross-section per grain for PAHs, VSGs and BGs. 
Fk,i{E) is the radiation field (in erg s^^ cm^^ eV^^), 
found from the transfer prescriptions described in Sect. 
3.3, and ni/(fc,j) the density in each cloud zone charac- 



terised by the index k ai id i (see Fig. 10). For the BGs we 
tak e the cross section of Desert et al. (1990 ) , for PAHs we 



use Verstraete fc Leger (1992 ) and Drainc fc Lee (1984 ) 



Table 1. Average power radiated per hydrogen atom by 
gas coolants and small dust grains towards L1721 





I 1 ricoTTrtiH 
V y Ubci V Kid 


' 1 ' Vi la/^Tdt 1 ci 1 
_L licUl t: LlUdl 


[C+] 158Mm 


1.2 


1 


[O"] 63Atm 


0.34 


0.32 


Ha 




0.27 


Total gas cooling 


1.81*' 


1.59 


IpAH + IvsG 


78 


54 


epE 


0.023 


0.029 



" In 10-2^ erg s'^ H"! 
'' Including the H2 contribution as estimated from the model 



for VSGs. In the ISRF, PAHs absorb 2.3 10"" Watt/C 
whereas graphitic VSGs only absorb 1.7 10"^"^ Watt/C. 
Formula (4) can then be inverted to infer the dust abun- 
dance (C/H in PAHs and VSGs) as a function of cloud 
radius. The resulting abundance profiles are presented in 
Figs. |ll|, |l^ and |l3|. The abundance of PAHs has to be mul- 
tiplied by ~ 5 in the [0.25 < Ay < 0.45]-layer to account 
for the observed Ipah/Ibg ratio. Similarly, the VSGs 
abundance must be multiplied by 2 in the [Ay > 0.75]- 
layer in order to reproduce the observed IvsgI Ibg emis- 
sion ratio. 

Including these abundance variations in the PDR mod- 
elling, we obtain the result displayed in Fig. |ri|. The cor- 
responding line intensities are higher: the C+ and 0° lines 
emission are increased by a factor about 1.5 and 3 respec- 
tively. The discrepancy at the center of the cloud between 
the observed and predicted line profiles may result from 
uncertainties in the radiation field, the radiation transfer, 
or the photoelectric efficiency and from our rough esti- 
mate of the small dust grain abundance profiles. In the 
present PDR modelling we have not attempted to include 
changes in the UV extinction curve associated with the 
small dust grain abundance variations because recent ob- 



servations (Boulanger et al., 1994) have questioned the 



contribution of small dust grains to the extinction curve 



adopted in the current dust models (Desert et al., 1990; 
pwck et al., 1997 , e.g. ). In any case, exploring a range 
of far-UV to visible extinction ratio (from 5 to 10 at 1000 
A) we find the cooling lines to vary by less than 50 %. We 
therefore conclude, that the observed intensity profile of 
the main gas cooling lines can only be reproduced if the 
abundance of PAHs is significantly increased. This result 
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illustrates the dominant role of the smallest dust grains 
in the gas heating. It also shows that the model photo- 
electric efficiency is close to its actual value. In Table 1, 
the values of the power radiated by gas coolants and small 
dust grains as observed and modelled are listed. 



4. Dust photoelectric efficiency 

Using our model estimate of the II2 cooling, we can now 
estimate the photoelectric efficiency by comparing the ob- 
served dust and gas emissions. A relationship is theoret- 
ically expected between the main cooling rates and the 
dust emission. Indeed, a dust grain exposed to a stellar 
radiation field dissipates the energy it absorbs mostly as 
IR emission. A small fraction, e, of the energy absorbed 
is channelled to the gas by the photoelectric effect. The 
heating efficiency of this latter process is defined as 

abs 

where Pabs is the total power absorbed by the dust 
grains and Ppe is the power given to the gas through the 
photoelectric effect. For a gas in thermal balance, the gas 
heating and cooling rates are equal. In the L1721 PDR, 
the gas heating is dominated by the photoelectric effect. 
Assuming thermal balance the photoelectric heating is 
thus roughly equal to the total cooling rate. As long as 
the opacity of the emerging cooling lines is small, we can 
replace PpE by the total gas cooling hue emissions /a. 
Moreover, as the power conveyed by the photoelectric 
effect is much smaller than the absorbed power (e ~ a few 
10-2), the latter power is well approximated by the dust 
IR emission. Finally, e can be observationally defined as: 



/a 

I dust 

where Idust is the total dust IR emission. 



(7) 




Fig. 14. Upper panel: intensities in the and CP lines 
measured by LWS, and cooling in the H2 rotational lines 
(see Sect. 3.3.4)- Lower panel: The total gas cooling, I\ — 
Iq+ +Iqo+I[{^ with 3a error bars (solid line with crosses) 
and as derived from the dust IR emission: (i) if^ — 0.03 x 
IpAH + 0.01 X IvsG + 0.001 X Ibg (solid line) and (ii) 
If^ = 0.001 X IpAH + 0.001 X IvsG + 0.02 x Ibg (dotted 
line). 



As the photoelectric heating rate is expected to be 
dominated by the small grain populations, it is interesting 



to determine the photoelectric efficiency for each grain 
population. Along the LWS cut, we define: 

^A = ^PAH IpAH + eySG IvSG + ^BG IbG (8) 

with epAH, ivsG and epG the photoelectric efficiencies of 
the PAHs, VSGs and BGs respectively and Ipah, Ivsg 
and Ibg the corresponding dust emission. The photoelec- 
tric efficiencies can then be found from a least-square fit 
of to I A- In this derivation, we assume a constant pho- 
toelectric efficiency for each dust component across the 
cloud: as we show below (see Fig. 15), this is a reason- 
able approximation. From our best fit, we find that the 
small grain populations have a high photoelectric effi- 
ciency, namely epAH — 3%, evsG — 1% and epG — 0.1% 
(see Fig. |l|). 

This combination reproduces well the observed gas 
emission in the outer parts of the cloud, but there is a 
significant deviation in the central region of the cloud. In 
fact, as discussed in Sect. 3.4, the gas cooling hues become 
optically thick towards the cloud center whereas the dust 
emission remains optically thin throughout the cloud: this 
is probably why clearly overestimates 1^ at r ^ 0.2 pc. 
We also show on Fig. |l^ the opposite case where the BGs 
dominate the photoelectric heating rate: the correspond- 
ing cooling emission I^^ is not a good representation of 
I A along the whole LWS cut. Our results clearly show that 
the dust photoelectric heating is dominated by the smaller 
grains, and not by the BGs. However, the spatial correla- 
tion between dust and gas emission does not allow us to 
tightly constrain the photoelectric efficiency of the vari- 
ous grain components because of the spherical geometry 
where outer and inner regions contribute to the emergent 
intensity at all sky positions. 




Fig. 15. Comparison between the theoretical efficiency 
£mod (solid line and see text) and our observational value 
e-obs = Ia/{Ipah + Ivsg) with Ik = Ic++Io'>+Ih2 (solid 
line with squares ) as a function of r and Ay . 

We can estimate the observational photoelectric effi- 
ciency profile of small grains across the cloud as eots = 
I a/ (I pah + Ivsg)- This efficiency is compared in Fig. [l^ 
to the theoretical curve, emod which has been derived from 
our PDR modelling with the dust abundance variations 
discussed in Sect. 3.5. For consistency with the observa- 
tions, we define Cmod as the ratio between the emergent 
gas cooling lines and the power absorbed by small grains 
as computed in Sect. 3.5. We see on Fig. |l^ that emod re- 
produces the main features of eobs across the cloud and 
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that the emod- values are within 30% of those of Eobs ■ The 
shape of the emod-profile can be understood as follows. 
While penetrating into the cloud (r >0.3 pc), the UV ra- 
diation is attenuated and the average grain charge (posi- 
tive) decreases; consequently, emod increases because of the 
weaker Coulomb barrier for the ejection of electrons from 
the grain (BT). Towards the center of the cloud (r <0.3 
pc), the photoelectric efficiency decreases. This trend is 
the combination of two effects. On the one hand, the opac- 
ity of the cloud center to UV photons is larger than for 
visible photons: since the photoelectric effect is produced 
mostly by UV photons, emod diminishes towards the cloud 
center. On the other hand, the emod-drop at r <0.3 pc also 
reflects the behaviour of the npAH /nysG ratio (see Fig. 
[l^ ). Comparing emod and eobs iu detail, we see that the 
efficiency values at the cloud boundary are quite compa- 
rable. Moving into the cloud, eobs is systematically lower 
than emod- Also, the ridge and trough of the emod-profile 
are somewhat shifted with respect to eobs- These discrep- 
ancies probably result from the simplified geometry and 
radiation transfer we adopted for the L1721 cloud, as well 
as the uncertainty in the spectral distribution of the excit- 
ing radiation field (see Sect. 3.2). In this latter case, a de- 
tailed modelling of the dust infrared emission is warranted 
in order to consistently derive the exciting radiation field 
and the dust abundance profiles in the L1721 cloud. 

5. Conclusion 

We study the gas thermal balance and the dust IR 
emission across LI 721 a nearby, isolated molecular cloud 
heated by the v Sco star and by the BA star associa- 
tion of Upper Scorpius. The coohng of the gas is obtained 
with ISO-LWS measurements of the major cooling lines, 
[C"*"] 158/zm and [0°] 63/j,to . The dust IR emission is 
traced with IRAS data. 

We model the gas cooling lines with an updated version 
of the photodissociation region model of Le Bourlot ct al 
1993| ). The input physical conditions in this modelling 



are the density profile and the exciting radiation field. 
The density profile of the L1721 cloud is deduced from vi- 
sual extinction data {Ay at the center is ~3) . The exciting 
radiation field is estimated to be x = 5 to 10 times the 
radiation field of the Solar Neighbourhood from the ratio 
of dust emission to gas column density. 
The spatial correlation between the gas cooling lines and 
the IR emission of small grains (PAHs and VSGs of radii 
< 100 A) confirms the theoretical expectation that the gas 
heat budget, in the less excited PDRs, is dominated by the 
photoelectric effect on small grains. The photoelectric ef- 
ficiency measured from the ratio of the total gas cooling 
to the dust IR emission is in the range 2 to 3% in good 



agreement with current theoretical models (Bakes & Tie- 



lens, 1994; Wcingartner & Draine, 2001) 



Moreover, we find that the gas cooling lines of the L1721 
cloud cannot be explained if the abundance of small grains 



is kept constant across the cloud, equal to the abundance 
of the diffuse ISM. Rather, as indicated by the IRAS band 
ratios, our PDR model results matches the gas cooling line 
observations when the abundance of PAHs is enhanced by 
a factor of 5 towards the cloud center. 
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